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Abstract 

The impact of free charges on the local pressure on a charged ferroelectric domain wall produced 
by an electric field has been analyzed. A general formula for the local pressure on a charged 
domain wall is derived considering full or partial compensation of bound polarization charges 
by free charges. It is shown that the compensation can lead to a very strong reduction of the 
pressure imposed on the wall from the electric field. In some cases this pressure can be governed 
by small nonlinear effects. It is concluded that the free charge compensation of bound polarization 
charges can lead to substantial reduction of the domain wall mobility even in the case when the 
mobility of free charge carriers is high. This mobility reduction gives rise to an additional imprint 
mechanism which may play essential role in switching properties of ferroelectric materials. The 
effect of the pressure reduction on the compensated charged domain walls is illustrated for the 
case of 180° ferroelectric domain walls and of 90° ferroelectric domain walls with the head-to- head 
configuration of the spontaneous polarization vectors. 

PACS numbers: 77.80.Dj 
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I. INTRODUCTION 



Understanding the dynamics of domain wall motion is essential for the explanation of 
many phenomena in ferroelectric materials. In many models for the evolution of domains in 
ferroelectrics, researchers deal with systems where charged ferroelectric domain walls appear. 
Classical examples to be mentioned are the work by Landauer^ on the nucleation of 180° 
domains and that by Miller and Weinreicb^ on the sidewise movement of 180° domain walls. 
In these models, the domain wall dynamics is analyzed by minimizing the thermodynamic 
potential consisting of the depolarizing field energy, the energy associated with the crystal 
lattice polarization, the energy of the domain wall, and the energy supplied to the system 
by the external electric source. It is evident that in all systems where charged domain 
walls appear there exists a possibility of compensation of bound charges by free charges, 
which affects the aforementioned energies through the reduction of the net charge at the 
domain wall. It was already pointed out by Landauer- that the process of bound charges 
compensation will not only reduce the depolarization energy, but it will also reduce the 
energy supplied to the system by external electric sources. Since the latter energy actually 
represents the driving force for domain wall motion, one can expect that the appearance of 
free charges in the system will result in a reduction of local pressure on the domain wall. This 
may have a serious impact on the domain nucleation and the sidewise domain wall motion. 
The problem of the impact of bound charge compensation on the domain pattern evolution 
seems to be of special practical interest for ferroelectric ceramics where the appearance of 
charged domain walls is readily expected whereas the high temperatures used to pole it can 
promote free-charge screening effects. 

The aforementioned issues have motivated the theoretical analysis presented below, where 
we will address the effect of bound charge compensation on the local pressure exerted on the 
walls by the electric field. We will cover the general situation as well as some cases of special 
interest. In Sec. [Ill we will illustrate the effect of charge compensation for the simplest 
situation. Section IIIII presents a derivation of the general expression for the local pressure 
imposed by the electric field to an arbitrarily compensated ferroelectric wall. In Sec. IIVI 
we will demonstrate the application of the general results to systems with low electric fields 
(Subsec. IIV Ap and to systems where nonlinearity of the dielectric response of the crystal 
lattice is essential (Subsec. IIVBII . 
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FIG. 1: Scheme of a ferroelectric capacitor of thickness h with a 180° ferroelectric domain wall 
tilted by an angle a with respect to its electro-neutral orientation. Bottom electrode is grounded 
and top electrode is at a constant electric voltage —U. Pluses indicate the positive bound charges 
(Tfi due to a discontinuity of the spontaneous polarization Pq at the domain wall, minuses indicate 
the free charges aj that fully compensate the bound charges. Symbol 5u stands for the virtual 
displacement of the domain wall and symbol p stands for the external mechanical pressure that 
should be applied to the wall to keep it at rest. 

II. SIMPLE MODEL WITH FULL CHARGE COMPENSATION 

At first we analyze the effect of free charges on the local pressure on a ferroelectric domain 
wall for the system shown in Fig. [H We consider the hard ferroelectric approximation where 
the electric displacement in it is presented as the sum of the constant spontaneous polar- 
ization ±Po with antiparallel orientation in neighboring domains and the linear dielectric 
response of the crystal lattice to the electric field, characterized by the permittivity e. 

We consider a planar 180° domain wall, which makes angle a with the vector of spon- 
taneous polarization. We address the situation where the tilt of the domain wall in the 
clockwise direction yields the appearance of positive bound charges of surface density 

= 2Pq sin a as it is shown in Fig. [H The bound charges are completely compensated 
by free charges of surface density o"/ = —cTb- The bottom electrode is grounded and the 
top electrode is connected to a source of constant voltage —U. The voltage applied to the 
electrodes produces an electric field Eq = U/h. 

In this study, we calculate the value of local pressure on the ferroelectric domain wall 
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FIG. 2: Detail of the intersection area of the domain wall and the top electrode indicating the charge 
variation on the top electrode during the virtual displacement of the domain wall. Dashed and solid 
lines show the original and displaced positions of the domain wall Sw, respectively. Symbol aAP 
stands for the surface charge density of bound charges due to spontaneous polarization reversal 
during the domain wall shift. Symbol crj stands for the surface charge density of free charges at 
the domain wall. Symbol 6Sap indicates the area where spontaneous polarization reversal occurs 
during the domain wall shift. Symbol 6Sf indicates the area of domain wall, which reaches the 
electrode during the domain wall shift. 

produced by the electric field using the principle of virtual displacements. We consider the 
wall at rest, when the force imposed on it by the electric field is balanced by the external 
mechanical pressure p. According to the principle of virtual displacements, the variation of 
the proper thermodynamic function 6G equals the work produced by the external mechanical 
pressure p during the virtual displacement, i.e., 6Wp = —SwP^u, where Sw is the area of the 
domain wall and 6u is the virtual displacement of the domain wall in the normal direction. 
One readily checks that for the considered model in the hard ferroelectric approximation, the 
thermodynamic function of the system consists of the electrostatic energy and the energy 
supplied to the system by external electric sources; it reads: 



where e is the permittivity of the ferroelectric and Q is the charge on the top electrode. 
Figure [2] shows in detail the intersection area of the domain wall and the top electrode. The 
considered virtual displacement is indicated by the line AC, where A indicates a point on 
the domain wall in the original position and C indicates the same point on the displaced 




(1) 



4 



domain wall. Since the bound charges on the oblique domain wall are considered to be fully 
compensated by free charges, the electric field in the ferroelectric capacitor is homogeneous 
and equal to Eq = U/h. This means that the variation of the thermodynamic function, 
6G, which is produced by the virtual displacement 6u of domain wall, is given only by the 
variation of charge on the top electrode, 6Q. Thus, we can write the following equation of 
equilibrium: 

U6Q = -SwpSu. (2) 

Considering the ferroelectric capacitor of length L in the direction perpendicular to the 
cross-section shown in Fig. [H the area of oblique domain wall equals 

Sw = Lh/ cos a. (3) 

The variation of charge on the top electrode 6Q includes the charge 6Q^p due to the spon- 
taneous polarization reversal, which is partially compensated by the free charge 6Qf that 
has arrived at the top electrode from the "annihilated" part of the wall during its vir- 
tual displacement 6u in the normal direction. Thus it is 6Q = 6Qap + SQf. The charge 
SQap is provided by the bound charge density ctap = 2Pq due to spontaneous polariza- 
tion reversal, which occurs on the area 6Sap on the top electrode indicated by the line 
BC in Fig. [21 Using trigonometric functions one can readily get the relation between 
the virtual displacement of the domain wall 6u and the area 6 Sap = L6u/ cos a. Thus 
the spontaneous polarization reversal produces the variation of charge on the top electrode 
SQap = —capSSap = —2PqL5u/ cos a. The charge 5Qf equals the free charge from the 
piece of the domain wall that "annihilates" at the electrode during the virtual displacement 
of the domain wall. This charge has a surface density of = —2Pq sin a and is distributed 
on the area 5Sf on the domain wall indicated by the line AB in Fig. [21 Using trigonometric 
functions, it is easy to express the relation between the virtual domain wall displacement 5u 
and the area 5Sf = —af6Sf = 2LPo'^'w(sina;)/(tana). Finally, one can present the variation 
of the charge on the top electrode in the form: 

6Q = -2PqL (l - sin^a) 5u/ cosa. (4) 

Substituting Eqs. © and into ([2|) we get the formula for the local mechanical pressure 
on an inclined 180° domain wall, which is fully compensated with free charges: 

p = 2Po^o(l-sin2a). (5) 
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FIG. 3: Special case of a fully compensated domain wall, which is parallel to the electrodes (i.e. 
a = 90°). In this case, shift 6u of the domain wall Sw in the normal direction does not produce any 
change of charge on the top electrode, 6Q = 0, resulting in zero local pressure on such a domain 
wall. 



It is seen that, for the domain wall that is normal to the electrode, i.e. a = 0, the above 
formula reduces to the well-known value of pressure 2PqE on an uncharged 180° domain 
wall. Figure [3] shows the opposite extreme case of a fully compensated domain wall, which 
is parallel to the electrodes (i.e. a = 90°). This case corresponds to a fully compensated 
head-to-head configuration of the spontaneous polarization. In this situation, according to 
Eq. the local mechanical pressure on the domain wall is zero. This can be readily 
understood. Indeed, in this case, the shift of domain wall 6u does not produce any change in 
the thermodynamic function G, since a displacement of the fully compensated wall should 
affect neither the electric field nor the charge on the top electrode. 



III. LOCAL PRESSURE ON A FERROELECTRIC DOMAIN WALL 

In this section we obtain a formula for the local pressure on the ferroelectric domain 
wall for the most general situation where (i) the electric field is inhomogeneous, (ii) wall 
configuration is arbitrary, (iii) charge compensation of the wall is not assumed to be full, 
and (iv) the analysis goes beyond the hard ferroelectric approximation. The only limitation 
of the presented analysis would be the neglect of the mechanical effects, i.e. we set the 
mechanical stress t^'^'^^ zero everywhere in the sample. Whenever such effects are essential 
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they may be incorporated in the framework. 

The desired formula can be obtained by apphcation of the virtual displacement method 
to a general configuration of electrodes and arbitrary compensated ferroelectric domain 
walls (see Appendix |A]). However, the calculations can be essentially simplified by using an 
advanced thermodynamic result derived in the classical textbook by Landau and Lifshitz^. 
According to this book, the volume force density fi in a dielectric can be presented in the 
form 

in term of the generalized stress tensor Tij which, in turn, can be written as 

= $ 6,, + r-^^^ + ^ {E,D, + E,D,) (7) 

with 

$ = $ - (1/2) 6oE^ - EiPi (8) 

where Pi and Di = Pi + e^Ei are the polarization and electrical displacement, respectively; 
$ is the thermodynamic function of the dielectric characterized by the differential = 
Ei dPi + T^"^"^^ duij where Uij is the strain tensor. 

Applying Eq. ([6]) to the two domains in the vicinity of a domain wall, one presents the 
pressure to the latter p in the form 

P = [[nj]] riiUj (9) 

where Ui is vector normal to the wall. The notation [[Z]] here and therein is used to denote 
the jump of the variable Z at the wall. Considering the mechanically free ferroelectric 
sample, i.e. r™^'^'^ = 0, and combining Eqs. ([7]) and one gets 

p=m] + [[E,D,]]n,n,. (10) 

If we further employ the expressions for the continuity of tangential components of the 
electric field at the domain wall [[-Et.i]] = [[E^ — (i?fcnfc)nj]] = 0, the relation for the normal 
component of electric displacement at the domain wall [[Di]] rii = af, and the algebraic 
identity 

[[E,D,]]=%[[D,]] + [[E,]]D,, (11) 



where Ei = {El'> + El'>)/2 is the average electric field at the opposite sides of the domain 
wall, we get the formula for the external mechanical pressure on the domain wall 

p = [[$ + (1/2) BoE' + P,E,]] - E, ([[P,]] - ajn,) . (12) 

Considering Eq. ([8]), formula ( fT2|) can be further simplified: 

P = - - I (i^f ^ + (Pf ) - Pf ) - . (13) 

Here upper indexes (1) and (2) specify to which domain the variable is referred to, the 
direction of the normal vector Ui is taken inside domain (2), and the pressure is considered 
as positive when acting from domain (2) to domain (1). 



IV. PRESSURE ON THE WALL IN SPECIAL CASES 

In this section, we present consequences of Eq. fll3p and demonstrate its applicability for 
two particular cases. At first we present further a simplification of formula fll3p for the case 
of low applied fields where the dielectric non-linearity of the ferroelectric can be neglected. 
Then we demonstrate the application of formula (JT3l) in systems where the non-linearity of 
the dielectric response of the crystal lattice is essential. 



A. Hard ferroelectric approximation 

For small electric fields we can use the "hard ferroelectric" approximation and express the 
polarization as a sum of the constant spontaneous polarization Po^j and the linear response 
of the polarization of the crystal lattice to the electric field 

P^ ~ Po,i + Xr.Ej, (14a) 

where Xij is the susceptibility of the ferroelectric. In this situation, function $(P) can be 
expressed in a form 

$(P) ^ ]pc^jE,Ej. (14b) 



After substitution of Eqs. (1141) into Eq. (1131) one readily gets 



V = -E. mA] - ^jn.) - \ [[Xd] E'^^'^Ef. (15) 
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FIG. 4: Normalized pressure on 180° ferroelectric domain wall versus the tilt of domain wall a with 
respect to vectors of spontaneous polarization, which is shown in the inset, for different degrees of 
compensation of bound charge ai, by free charge af, A = —af/ab- 

In the case of absence of free charges in the system, Eq. ( fTSl) was obtained earlier-i^. 

Two numerical examples of the effect of free charges on the local pressure on a charged 
ferroelectric domain wall are shown in Figs. H] and [5l In Fig. H] we consider the case of 
a 180° ferroelectric domain wall where the vector of average electric field Ei is considered 
to be always parallel to the vectors of spontaneous polarization P^ i and the domain wall 
is considered to be tilted by an angle a with respect to the directions of spontaneous po- 
larization. The tilt of the domain wall results in the appearance of bound charge on the 
domain wall, cxb, and the local pressure of external sources on such a wall equals 2Po,i-E'i 
and its direction is shown in the inset of FigJH Then we consider that the bound charge is 
partially compensated by free charges af. The degree of partial compensation is expressed 
by the quantity A, so that Cf = — Actj, where < A < 1. It is seen that free charges reduces 
the value of pressure p, but the important point is that, for the same degree of screening A, 
the value of p varies depending on the domain wall tilt. Zero pressure corresponds to the 
situation of full compensation of the domain wall, which is perpendicular to the vectors of 
spontaneous polarization. 

The effect of free charge compensation on the pressure on a charged 90° domain wall is 
shown in Fig. O However, the contribution due to dielectric anisotropy is neglected (i.e. 
[[Xij]] = 0), since it produces a field redistribution^, which can be hardly evaluated in the 
general case. Since this type of domain wall is ferroelastic it is reasonable to consider that 
the angle between the vectors of spontaneous polarization and the domain wall is fixed as 
it is indicated in the inset of Fig. [5l The domain wall charge is equal to at, = V2Po and 
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FIG. 5: Normalized pressure on a 90° ferroelectric domain wall versus the angle, a, between the 
normal to the wall and the vector of electric field, which is shown in the inset, for different degrees 
of compensation of bound charge ai, by free charges aj; A = —Of/a},. Note that for complete 
screening of bound charges, i.e. A = 1, the local pressure on the domain wall is zero for any 
relative orientation of the domain wall and the electric field. Here the contribution due to the 
dielectric anisotropy is neglected, i.e. [[xij]] = 

the local pressure is a function of the tilt of the domain wall with respect to the average 
electric field Ei. It is seen that, in this case, the presence of free charges leads to an effective 
decrease of the value of spontaneous polarization. The physical reason for this is that in this 
case the vector [[-Po,i]] is always parallel to the normal vector of the domain wall n^. It is 
noticeable that, in the case of the complete bound charge compensation, the pressure on the 
domain wall is zero regardless the orientation of the average electric field and the domain 
wall. 

B. Nonlinear effects 

Formula fll3p makes it possible to analyze the effects of the dielectric nonlinearity of the 
crystal lattice on the local pressure p. As an example we consider a fully compensated 180° 
domain wall, which is perpendicular to the vectors of spontaneous polarization (see Fig. [2]) . 
In section |TT] we have shown that, in the hard ferroelectric approximation, the local pressure 
on such a wall is zero. However, analysis beyond this approximation shows that this force 
is not actually zero. Let us find this force for the case of a uniaxial ferroelectric with the 
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second order phase transition. In this case, the free energy of the ferroelectric reads: 

$(P) = $o + ^«P^ + ^/?P', (16) 

where, in the ferroelectric phase, a < and f3 > 0. Discontinuous change of polarization 
states P*-^-* and P*-^-* at the domain wall and the presence of free charges o"/ produce the 
electric field AE according to the equation: 

p(2) _ p(i) _ = 2eoAE. (17a) 

This electric field is of antiparallel orientation at the opposite sides of the domain wall and 
it is superposed on the uniform applied electric field Eq, so that the equation of state in the 
adjacent domains are of a form: 

5$(p(i))/5P = Eo + AE, (17b) 
d^{P^'^'^)/dP = Eo- AE. (17c) 

The polarization P^^^ and P*^^) at the opposite sides of the domain wall is given by the 
solution of Eqs. (fT7|) and, in the case of fully compensated domain wall, i.e., o"/ = —2Pq, it 
can be expanded in a Taylor series with respect to the applied electric field Eq: 

^'^'--^«H-.^.„.^ji|^i.„^ + ..., (18b) 



where Pq = y'—P/a is the spontaneous polarization and x = ~V(2ct) is the dielectric 
susceptibility of the crystal lattice. After substitution of Eqs. ffTSl) into Eq. f|T6l) . we can 
readily obtain the expansion of function $(P) at the opposite sides of the domain wall in 
terms of the applied electric field: 

When one considers that the average field in the ferroelectric, E, is dominated by the applied 
electric field Eq, i.e. E = Eq, the resulting formula for the pressure on the domain wall is 
according to Eq. (fTHIl equal to: 

P ~ - ixVPo) El (20) 
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FIG. 6: The applied field dependence of the local pressure on a fully compensated head-to-head 
180° ferroelectric domain wall. The pressure is normalized to —2PqEc, that equals the value of 
pressure expected to be applied to an uncompensated domain wall at the thermodynamic coercive 
field Ec in the hard ferroelectric approximation. The field is normalized to the thermodynamic 
coercive field Ec- The exact values (solid line) obtained by numerical solving the nonlinear Eqs. 
p7|) are compared with the approximative values (dashed line) calculated using the formula (j20|) . 

which means that the pressure on the domain wall is oriented in the opposite direction than 
the "expected" 2PqEq contribution. 

The physical interpretation of the above formula is as follows. The application of a 
high electric field to the ferroelectric sample with the considered configuration of a fully 
compensated charged 180° domain wall is associated with two effects: first, it produces the 
difference in volume energy density of the crystal lattice at the opposite sides of the domain 
wall due to dielectric nonlinearity of the crystal lattice, and second, it produces an additional 
charging of the domain wall due to a different polarization response at opposite sides of the 
domain wall so that the wall becomes negatively charged. The pressure associated with the 
electrostatic force acting on this charge is directed against that associated with the energy 
difference. Finally, the resulting formula corresponds to the sum of these effects. 

Figure [6] shows the applied electric field dependence of the local pressure on a fully 
compensated domain wall, which is parallel to the electrodes (see Fig. [3]). The dashed 
line shows the approximative result given by formula ( 120|) . The solid line shows the exact 
result given by the numeric solution of the system of nonlinear equations f|T7|) . where the 
function $(P) is given by Eq. fll6p . In Fig. [HI numerical values typical for barium titanate 
were used: Pq = 0.26 Cm~^ and x = 188 ^o- It is seen that Eq. fl2U|) provides a very good 
approximation, which differs only by about 6% from the exact value at the thermodynamic 
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coercive field Ec = Pq/ (3-\/3x). 



V. ADDITIONAL IMPRINT MECHANISM 

The results obtained in the previous sections enable us to draw a conclusion on the 
imprint behavior of polydomain ferroelectrics. 

Imprint is a property of a ferroelectric capacitor of exhibiting a higher switching voltage 
when the switching starts from a state in which the capacitor has been kept for a long time 
or subjected to light illumination (see Ref.- for comprehensive discussion). A popular and 
widely referred to imprint mechanism, which has been offered by the Sandia Laboratory 
group®'^, is related to free-carrier compensation of charged domain walls. Our analysis 
presented in this paper suggests that this compensation should lead to one more imprint 
mechanism. 

The key issue of the Sandia mechanism, which we will refer to as "trapping imprint", 
is that, if the free carriers compensating the bound charge of a wall are trapped, then 
the mobility of this wall is deteriorated since now it is to drag with it the weakly mobile 
trapped carries. Since, in typically low-conductive ferroelectrics, an appreciable time or/and 
illumination is needed for the free-carrier compensation of the bound charge of a wall, 
this effect will lead to the imprint. In the context of the analysis presented above its is 
clear that the aforementioned charge compensation will also lead to a reduction of the 
"thermodynamic" (pondermotoric) force applied to the wall in the switching field. This 
reduction, similarly to the above case, can manifest itself in an imprint mechanism, which 
we will refer to as "force- reduction" imprint. One should stress that the performance of these 
mechanisms can be appreciably different. For instance, in the case of trapping imprint, if 
the switching field is applied long enough to let the trapped charges drift with the wall, 
the imprint effect becomes weaker. Its efficiency is clearly controlled by the mobility of the 
trapped free carriers; in the limit of high mobility the effect vanishes. On the other hand, 
the force-reduction imprint is clearly insensitive to the mobility of the compensating free 
chargers. In contrast to the trapping imprint, one cannot cope with this mechanism by 
using longer switching pulses. All in all, one sees a possibility to distinguish these imprint 
mechanisms in real ferroelectric systems. 
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VI. CONCLUSIONS 



The impact of free charges on the local pressure on a charged ferroelectric domain wall 
produced by an electric field has been analyzed. A general formula for the local pressure 
on a charged domain wall is derived considering full or partial compensation of the bound 
polarization charge with the free charge. It is shown that the compensation can lead to 
a very strong reduction of the pressure imposed on the wall by the electric field. In some 
cases this pressure can be governed by small nonlinear effects. It is concluded that the 
free charge compensation of the bound polarization charge can lead to substantial reduction 
of the domain wall mobility even in the case of the high mobility of the free charge. This 
mobility reduction gives to an additional imprint mechanism which may play essential role in 
switching properties of ferroelectric materials. Using this mechanism, it is possible to explain 
frequent observations of very stable domain patterns with the head-to-head configuration 
of the vectors of spontaneous polarization^^. The results obtained can be used in modeling 
the poling of ferroelectric ceramics, domain nucleation, sidewise domain wall movement and 
spontaneous polarization switching in imprinted, leaky ferroelectric samples, and samples 
treated at elevated temperatures. 
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APPENDIX A: LOCAL PRESSURE ON A FERROELECTRIC DOMAIN WALL 
DERIVED FROM THE PRINCIPLE OF VIRTUAL DISPLACEMENTS 

In this Appendix we offer an alternative derivation of the local pressure on ferroelectric 
domain wall than it is presented in Sec. Illli Here, we use again the principle of virtual 
displacements as we do in Sec. [Tll We consider a system shown in Fig. [7] where inside a 
ferroelectric material with polarization Pj^^^ there exists a domain of the material with a 
different polarization state Pj , i.e. that there is a discontinuity of the polarization across 
the interface Sw This interface splits the ferroelectric into two domains, which have volumes 
V^-^^ and V^'^\ respectively. Inside each domain there is a conductor, which carries charge 
Q^^^ and Qe\ respectively, and has electric potential ip^^^ and ip^^\ respectively. We consider 
that bound charges due to discontinuous change of polarization at the domain wall Sw are 
partially compensated by free charges of surface density o"/. The charges on conductors and 
on the domain wall produce electric fields E^^^ and E^"^^ within each domain. Symbols y?'-^^ 
and (f^"^^ stand for electric potentials within each domain. 

In order to obtain a general formula for the external mechanical pressure p on the domain 
wall, we follow the principle of virtual displacements presented in Section [Tll However, we 
adopt the thermodynamic function G in a more general form: 



The first term in the above formula represents the energy associated with the polarization 
of crystal lattice of ferroelectric, $(-P), and the electric field energy, {l/2)eoE^, which are 




(Al) 



15 



FIG. 7: Configuration used to calculate the local pressure imposed by the electric field to a wall. 
Domain wall Sw separates the ferroelectric into two domains of volumes V^^^ and V^'^\ There 
are two conductors, which carry charges and and have electric potentials ip^^^ and </Jg^ 
within each domain. The bulk quantities within each domain are the polarization pj^^^ and Pj^'^\ 
electric field E^^^ and E^'^^ and electric potential ip'^^^ and ip^'^\ Symbol 5u stands for the virtual 
displacement of the domain wall and symbol p stands for the external mechanical pressure that 
should be applied to the wall to keep it at rest. 

integrated over the volume V of ferroelectric. The second term in Eq. (lAip represents the 
subtracted work of the electric sources. It should be also noticed that the thermodynamic 
function G given by Eq. ( lAlj) is used just for convenience; the result (equation of state) is 
independent of the choice of the potential as always. 

In what follows it is convenient to transform the work of electric sources into volume 
integrals: 

y,«Q« = -y,g)/ Df\aS=! Ef'^DfUv+j ^^'^Df\aS, (A2a) 

i4'> JvW Jsw 

^'S Q'S = -^f I Df\aS = I Ef^DfUV-l ^^^^D^n^dS, (A2b) 

Js'-^^ Jv(^) Jsw 

where D^^^ and D^^^ are the vectors of electric displacement in domains V^^^ and V^'^^ 
respectively. In the derivation of the above equations, we consider the absence of free 
charges, i.e. divDj = 0, inside domains V^^^ and V^"^^ and vanishing the surface integral 
over the domain V^^^^ at infinity. Orientation of the normal vectors rij is indicated in Fig. [71 
If we further employ the expressions for electric displacement Di = Pi + eoEi, the continuity 

16 



of the electric potential at the domain wall [[</)]] = 0, the continuity of the normal component 
of electric displacement at the domain wall [[-Dj]] rti = af, and the algebraic identity used 
in Eq. ( fTTl) . the thermodynamic function G can be expressed in the form: 

^(1)77.(1) 77.(1)0(1) 



G = [ <^^\P^)-leoE^^E^-El^p[- 



dV + 
dV + 



(A3) 



+ / (pafdS. 

'Sw 

For the sake of presentation, it will be useful to define functions 

$B = $(P) - ^eoE.E, - EiPi, 
= ^o-f 

and to write the function G in the form 
G= [ ^BdV+ [ ^sdS, 

Jv Jsw 



(A4a) 
(A4b) 

(A5) 



where the volume integral is taken over the volume V = V^^^ + V^^\ 

Now our task is to express the variation of thermodynamic function 6G, which is produced 
by the virtual displacement of domain wall: 



5G 



dV 



V 



S<^s dS - 

Sw J Sw 



B 



dxi 



■Uk } 5u dS, 



(A6) 



where the first two terms represent the variations 6^b and S^s due to a change in bulk 
quantities during the virtual displacement of domain wall 6u. The last term represents the 
contribution to the variation 6G due to the volume change of domains produced by the 
virtual displacement of domain wall 6u and due to the change of the position of the domain 
wall. At equilibrium, the variation 6G equals the work 6Wp produced by external mechanical 
pressure p during the virtual displacement of the domain wall 



SWr, 



p6u dS. 



(A7) 



Sw 



Employing the principle of virtual displacements 6G = 6Wp , it can be readily shown that the 
first two terms in Eq. ( ]A6[) vanish because of the equations of state in the bulk ferroelectric 
(9$ 



OR 



Ei, 



(A8) 
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Gauss' law for electric displacement dDi/dxi = 0, and the continuity of the normal compo- 
nent of electric displacement at the domain wall [[-Dj]] rii = af. Finally, the last term in Eq. 
(IA6I) yield the formula for the external mechanical pressure on the domain wall 

-ElafUi. (A9) 

Taking into account the continuity of tangential components of the electric field at the 
domain wall, [[-Et.i]] = [[Ei — {Eknk)ni]] = 0, continuity of the electrostatic potential at the 
domain wall [[ip]] = 0, and the algebraic identity (|TT1) . we thus arrive at the result announced 
above, Eq. (fT3|l . i.e.: 

p = [[^p)]]-%m]]-^fn^). (Aio) 



$(P) - -eoEiEi - EiPi 
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